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^ , Abstract 

o 

P^ , We study permutation type solutions to n-simplex equations, that is, solutions 

^-^ ■ whose matrix form can be written as i?f^'"/" = Y\lt-^ ^''"n with some n x n 

^^ I matrix A and vector B, both over Iajj. With this ansatz the iD"("+^) equations of 

,__^ i the n-simplex equation reduce to an [in(n + 1) + 1] x [in(n + 1) + 1] matrix equation 



> 

X 



2 "v" I -^; I ^J L2' 

over Z/). We have completely analyzed the 2-, 3- and 4-simplex equations in the 
generic D case. The solutions show interesting patterns that seem to continue to 
still higher simplex equations. 

1 Introduction 

The Yang-Baxter equation (YBE, or 2-siniplex equation) is the fundamental equation of 
solvable models in (l+l)-dimensions. For lattice models it guarantees the commutativ- 
ity of the transfer matrix, and for particle scattering it implies solvability through the 
factorization of the scattering matrix [|l],|2|. Therefore, in order to construct interesting 
solvable models one needs interesting solutions. For this reason the YBE has been studied 
extensively and indeed many solutions are known [0,|^, especially in the two-state case |^. 
When one tries to generalize these solvable models to (2+l)-dimensions, either by con- 
sidering 3-dimensional lattices or the scattering of straight strings, one obtains Zamolod- 
chikov's tetrahedron equation (3-simplex equation) as the fundamental equation 
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whose solutions are needed for further development. Unfortunately only a few solutions 
are known for this equation ||^-|^ and when one proceeds to still higher dimensions and 
to the corresponding higher simplex equations very little is known. 

The difficulties associated with these equations come mainly from sheer numbers, the 
D-state n-simplex equation is actually a set of £)"("+^) equations on D^^ variables (in the 
non-constant case {n + 1)D^" variables). Because of this one is forced to make rather 
restrictive ansatze in order to obtain any solutions at all. One method is to take some 
definite high level structure (Lie algebra, chiral Potts) coming from somewhere else and 
apply it to the present situation. Our approach is complementary to this, the ansatz given 
below is defined in rather simple terms and we will then determine all solutions within 
this class. 

Let us recall the standard setup for the n-simplex equations. As usual we assume that 
we have linear operators TZ which for the n-simplex case are assumed to act on a product 
of n identical vector spaces V, i.e., TZ : V^®" -^ y^n^ j^^^ ^^ ^^ ^^le D basis vectors of V. 
Since we want to do algebra with the indices of the basis vectors it would be nice if the 
indexing formed a finite field. If D is prime this is possible with Z/j, integers modulo D, 
which is what we consider in this paper. However, some aspects of the following derivation 
works even if the indices just form a ring, for example with Z4. 

To the operator TZ we associate a numerical matrix with n pairs of indices by 

n (e,, ® ■ ■ ■ ® e,J = Ri^-::i: {e,, ® ■ ■ ■ ® e,J. (1) 

(Here and elsewhere in this paper summation over repeated indices is assumed.) The n- 
simplex equation itself is defined on V"®["-("+^)/^l, and the linear operators operate trivially 
in all but the n spaces indicated by the subscripts, e.g.,7li2 (cj^ e^j ® 623) = RlH^ (cj^ ® 
ej2 <S) Cjg), or in the general case with Ka G {1, . . . A^}, A^ = jn{n + 1), 

A^ 

fry \jl--JN _ Ty}Ki---JKn TT cjk fr)\ 

[RKi...K„)i,...iN - -tiij^^...iK„ 11 ^ik- l^J 

k^KcWa 

In this paper we consider the first few constant simplex equations, those given by the 
2-simplex or vertex Yang-Baxter equation 

7^l27^l37^23 = n2^nr^n^2, (3) 

the 3-simplex or tetrahedron equation 

'^123'^145^246^356 = '^356^246^145^123, (4) 



and the 4-simplex equation 

'7^1234'^1567'^2589'^3680'^4790 = '^479o'^368o'^2589'^1567'^1234- (5) 

In terms of the multi-indexed matrices defined in (H) the above operator equations imply, 
respectively, 

73^2^3 j^kih Tihh _ T3kik2 phka phh /fi\ 

^J2J3 ^jik-i -^fcifca ~ ^jij2 -^fciia -^fcafcs' ^^) 

Tjkzk^kd r)k2k4,lQ p^iU's jDhhh D^ifc2fc3 jDhk4k5 rthUke jDhhle ('7\ 

^J3J5J6 ^J2J4k6 ■'^jik4k5 ■'^kik2k3 — ■'^jiJ2J3 ^k-ijij^ ■'^k2k4,j,i-'^kzk^k,^i \' ) 

r)k4krkgko pfcafcefcgZo TDk2kc,lglg Tjkil^lQlr p'l'2'3k __ 
J4J7J9J0 jajejgko J2Jsksk9 -""Jlfesfcefc? -^kik2k3k4 ~ 

■nk\k2kzk4 r^lik^kfiky jDhhkakg -rihl&hko TjUhhlo (Q\ 

In addition to the above some other equations have appeared in the literature, e.g., the 



Frenkel-Moore equation [IC]. For a general formulation of the various type of equations. 



see 11 . 



2 Formulation with the permutation ansatz 

In this paper we consider only permutation type operators, that is those which transform 
one product of basis vectors into another simple product. In the matrix form this means 
that there is precisely one nonzero (= 1) entry in each column and row. For the In- 
state n-simplex equation there are {D"')\ different matrices to consider, and a brute force 
check of them is out of question except for D = 2, n = 2 which contains 24 permutation 
matrices (the next cases (2^)! = 40320 and (3^)! = 362880 might still be possible). We 
will therefore make the further assumption that the dependence between the basis vectors 
is linear, that is, 

7^ (cji (g)---0eij = eA<^i^+B^ ® ■ ■ ■ ® eA-i,+ij„, (9) 

(where the summation over a runs from 1 to n) for some nonsingular n x n matrix A and 
ra- vector B, both having entries from Z/). In terms of the i?- matrix this means that 

The main advantage of this ansatz is that the problem of solving the D-state n-simplex 
equation can be reduced to handling ordinary matrices over Zd, as will be shown below. 



This simplifies the problem considerably. Furthermore, although possible applications 
normally imply further conditions on the solutions, permutation matrices are such fun- 
damental objects that there is a good change they are acceptable in most cases, and we 
believe that the ansatz is not an unnatural starting point. 

In order to write the n-simplex equations in terms of A and B let us further define (in 
analogue with (||)) 

. A^, iii = Ka, j = Kp for some a, /3, 
[AK,...Kji = \ . ^, . (11) 

0^, otherwise, 

Bi, iii = Ko,, for some a, 
{BKi...Kji = < ^, . (l^J 

[ U, otherwise, 

so that 

N 

11=1 

where now the u summation runs from 1 to A^. 

In the homogeneous case, that is with B = 0, the above correspondence between R 
and A means that the n-simplex equation with ansatz (|1^) becomes an A^ x A^ matrix 
equation over Z^,. For example, the 2-simplex equation becomes 

(A2)f (Ai3)r(^23)L. = {A,s)'liA,,)^iA,,)l, (14) 

where A^l are 3x3 matrices with entries from Z^) as given in (|lID (for the explicit form 
see (0)). 

In the non-homogeneous case with B y^ a matrix formulation can also be obtained, 
if we add a fictitious index space and write 

n 
TDJl-jn _ dl . . _dv. do _ TJ d^ /I C-N 

-"■«!. ..i„ — "A«J„+Biio "A^io^+Bnio^io ~ 11 "A-^i^- ^ '-^^ 

fj,=0 ^ 

When this is immersed in the larger spaces we write the new index as the last one and 
then for the n-simplex case we get the [\n{n -|- 1) -|- 1] x [|n(n + 1) -|- 1] matrix 



A 



Kt...K„ 



Aki...K„ Bk^_k„ 

1 



(16) 
[We use square brackets when writing out these index matrices.] For example the 2- 



simplex equation becomes 



Al2 


Bl2 




"Ai3 




^is" 




^23 




-B23 







1 




1 




1 






^23 -^23 




^13 -Bl3 




A12 


B12 









1 









1 







1 



(17) 



and expanding this yields ( ]T^ ) for A and 

^12^13-823 + ^12-Bl3 + B12 



^23^13-Bi2 + ^23-Bl3 + B 



-"23; 



for B. The higher simplex equations have equally simple matrix form. In fact, formally 
the equations now look exactly as in (^§) with A instead of i?, but the interpretation is 
different: for A we have ordinary matrix products. 

3 Symmetries 

Before starting to solve the equations it is necessary to discuss their symmetries. For 
one thing, we only want to list the basic solutions from which the others are obtained by 
the allowed transformations. It is well known 0] that the ra-simplex equations are form 
invariant under discrete transformations of index transposition and index reversal. Now 
we should also see if these transformations preserve the linear permutation structure and 
what they imply on A and B. 



3.1 Index transposition of R 



If n^-t 






is also a 



is a solution of the A^-simplex equation, then {IR)ll'''l2 
solution. This is easy to see from the structure of the equation. 

Let us now see what the above symmetry implies for the index matrix A. From the 
definitions above it follows that 



iiR)tZ 



n ^'kjp+B^ 



w^i 



a=l 



{IA)iii3+{IB)^- 



(19) 



and by comparing the two expressions we find that if i? = 5{A^ B) is a solution, then 
[IR] = 6{A-\ -A-^B) is a solution, that is, (M) = A-\ (IB) = -A'^B. 

It is easy to see that this is also an invariance of the A equation. As a matrix equation 
it is clearly invariant under matrix inversion (which furthermore does not change the 
location of the inserted pieces of the unit matrix) and 



A Bl-1 
. 1 J 







-A~'B 
1 



3.2 Index reversal of R 

It is also easy to see that if RJlZi^^ is a solution then {CR)ll"'j'^ := Ri^Zll is a solution. 
We have 

n n 

(CHtt = n *^:i:..,«„,_. = n ^j;,,..,„c«.. (20) 

a=l a=l 

and the comparison yields (CA)^ = A^+iI^, {CB)a = i?„+i-a, that is, reflection across 
the center of the matrix or vector. 

Since the A-equation is a matrix equation it is invariant under any permutation of the 
set over which the summation is taken: if A^ solves the equation, then {(tA)^ := A'^.l 
where a is any permutation operator, is also a solution. However, we have to keep intact 
the structure of inserted parts of the unit matrix in the various terms, and then it appears 
that only the above reversal is possible. 

3.3 Transposition of A 

From the point of view of A the matrix equations have one more discrete symmetry: they 
are invariant also under transposition. However, this does not seem to correspond to any 
obvious invariance of R. In the following it will turn out that often the transposition 
of a solution A is also obtained by the central reflection (accompanied with parameter 
changes). However, this is not always true, and when it is not, it turns out that often the 
accompanying B will also be different. 

This is a rather interesting result from the point of view of studying the structure of the 
equations. Normally imposing ansatze on the solutions restrict the symmetries, because 
the symmetries of the equation may not be symmetries of the ansatze. In the present 
case this happens with the continuous transformation below. However, the opposite can 
happen as well: in the present case the ansatz leads to a new formulation which has its 
own obvious symmetries, and some of these do not seem to have any counterpart at the 
original level. 

3.4 Gauge transformations 

The gauge transformation Rki...k^ -^ {QR)ki...k„ = Qk\ ■ --QkI Rki...k„ Qki ■ --Qk^, is 
also an invariance of the n-simplex equations. Now that R is made out of delta-functions 
the transformation matrix Q must also be of that form, i.e.. 



(If D is not prime m ^ is not always defined.) A simple calculation shows that if R 
5{A, B) is a solution, then QR := 5{A, QB), where 



iQB)^:=uB^ + il-J2Al] 



(21) 



is also a solution. Thus only B can change, and we can in fact put one S^ = 0, if 
J^-yA"^ 7^ 1. Later we will find that for many solutions the inhomogeneous B part is such 
that it can be completely eliminated by this gauge transformation. 

In order to understand this as an invariance of the equations (|T^,0) we note first that 



(p!^) can be written as 



^12^13(1 - A23) + Au{l - An) + (1 - Au) = ^23^13(1 - Au) + ^23(1 - A13) + (1 - A23). 

If we now sum over the rightmost index of this equation and take its linear combination 
with equation (|18D we get (0) for (QB). 

Matrix equations are invariant under a much larger group of similarity transformations: 
A -^ O^^AO, but now that we have to preserve the structure of having inserted pieces 
of the unit matrix in Axj^...k„ these similarity transformation are allowed only with the 



matrix O 



xl yl 
1 



corresponding to the above. 



4 Results for the 2-simplex equation 

The details for the Yang-Baxter or 2-simplex case are as follows. In the homogeneous 
case we write 



pJiJ2 _ «i rj2 

ili2 ^aii+bi2 cii+di2'> 



(22) 



so that ^ = [c d] and then the 2-simplex equation becomes 



a b 




c d 




_ 1 _ 





a b 




1 




c d 





100' 




"100' 




a b 


= 


a b 




c d 




c d 





a b 




1 




c d 





a 


b 


" 


c 


d 











1 



(23) 

(In |12| similar matrices, but with the entries also being matrices, are used to define a 
dynamical system.) This yields 5 equations, 



abc = 0, bed = 0, bc{b — c) = 0, b{ad + b — 1) = 0, c{ad + c — 1) = 0, 



(24) 



whose solutions are discussed below. 



4.1 D is prime 



Recall that we are working with integers modulo D, the number of states. If D is prime 
there are no divisors of zero, and we can solve the equations with conventional rules of 
algebra. It is easy to show that in this case there are precisely four nonsingular solutions: 



A 



(1) 



a 
d 



A, 



(2) 



a 1 — ad 
d 



A, 



(2r) 



a 
1 — ad d 



A 



(3) 



1 

1 



(25) 



Whether these are really different depends on D. Note for example, that for D = 2 cases 
(2) and (2r) reduce to case (1); the same holds with a = d = 2 when D = 3. Solution 
(2r) is obtained from (2) by central reflection, and there is actually no need to mention it 
separately. 

When the inhomogeneous part B = [x,yY is included we have to solve ([T8|), which 
amounts to 

b{x + ay) = 0, c{y + dx) = 0, x{c + d — be — 1) = y{a + b — be — 1). (26) 



The solutions then split further and we get 



[A\B] 



[A\B 



(la) 



(2) 



1 
1 



[AlBf'^ 



a 1 — ad 
d 



—az 

z 



[A\B] 



a 
d 

(3) _ 



(a- l)z 
{d-l)z 



1 

1 



However, we have not yet used the gauge freedom. For [A|-B]2 " and [A|i?]2 the row 
sums of A are = 1, and thus according to ( pT]) we cannot change the inhomogeneous part, 
except by an overall multiplication. For [A|-B]2 the gauged inhomogeneous part will 
turn out to be 



{QB)i 

{QB)2 



(a — \){uz — v) 
{d — l){uz — v) 

(2) 



and by choosing v = uz we get {QB)i = 0. For [v4|i?]2 we get similarly 



{QB)i 

{QB)2 



-a{uz + f (1 — d)) 
uz + v{l — d) 

Now if (i 7^ 1 we can again transform to {QB)i = 0, but if rf = 1 only scaling is possible. 
Thus the final form of the solutions of the 2-simplex case is 



[A\Bf''^ 



1 


X 


, [A\Br^ = 


a 


' 


, imr = 


' 1 


" 


1 


y . 




d 







1 






a 1 — a 


—az 


, [m'? = 


a 1 — ad 


n 










1 


/C 




d 


U 



[A\B]^t^ 



up to the allowed transformations. 

4.2 ^ = 2,3 

For D = 2 the above yields basically two solutions for the 2-siniplex equation, R is either 
the unit matrix (with possible inhomogeneities) or the permutation matrix P, 



[A\B] 



' 1 


X 


or 


" 1 


' 


1 


y . 




1 






(27) 



where x, y G Z2. The same five solutions are obtained by a brute force search without the 
linearity assumption. 

The results (^^ work for any D, but already ioi D = 3 we get other homogeneous 
solutions, including triangular ones: 



1 
2 



") 


"20' 

2 


5 


"12" 
2 


5 


"22" 
1 


7 



2 2 
2 



and their reflections. 



4.3 D = 4 

The situation is quite different ii D = 4, because of divisors of zero: 2-2 = (mod 4). 
In addition to the above generic solutions, we get new base solutions 



"12' 

1 


; 


"12' 
2 1 


' 


"12' 
2 3 


) 


"32' 

3 


5 


"32' 
2 1 


' 


"32" 
2 3 



and their reflections. When the inhomogeneous parts are added we get 



" 1 


2 









1 


2 


y 


+ 2x + gx 




" 1 2 


C 


) 


1 


2x 


5 


2 1 


y + gx 


5 


2 3 


2x 




" 3 2 


2x ' 




" 3 2 


2x ' 




" 3 2 


y + 2x 









3 


C 


) 


5 


2 


1 





5 


2 


3 


1/ 







In the second solution there is no obvious way to fix the gauge parameter g and it has 
been left open. 

Thus for 4-state models there seem to be are additional symmetries and solutions, and 
perhaps this case needs more detailed studies. 



5 Results for the 3-simplex equation 

Higher simplex equations have many reductions to lower simplex equations, and it is not 

necessary to repeat them. For example, any solutions of the 2-simplex equation generates 

a solution of the 3-simplex equation by Rijk = RtjSk or SiRjk- These solutions (and those 

with det A = 0) will not be included in the following list and the solution below are 

genuine 3-simplex solutions. Note also that Rijk = Rik^j (with 6 on the central index) 

is not automatically a solution, in particular the permutation matrix RiH^lj = Sj'^Sj^SJl 

does not solve the tetrahedron equation. 

In order to solve the tetrahedron equation under the present ansatz we first consider 

the homogeneous part. The equation to solve is just like (^ with TZ replaced with A. 

When the matrix 

a b c 



A 



X y z 
u V w 



is inserted into the 6x6 matrix Akj^K2K3 the six different ways indicated in (^) and we 
compute the corresponding matrix product we find 29 equations: 

abx = 0, bxy = 0, vyz = 0, vwz = 0, 

bx{b — x) = 0, vz{v — 2;) = 0, y{bu — cv) = 0, y{—cx + uz) = 0, 

b{ay + b — 1) = 0, x{ay + x — 1) = 0, z{wy + z — 1) = 0, v{wy + t> — 1) = 0, 

abuz + acx + bcu = 0, bvxz + cvy + cxy = 0, buwz + cuz + cvw = 0, 

abu + acvx + cux = 0, buy + bvxz + uyz = 0, cuv + cvwx + uwz = 0, 

abwz + acz + bcw + c^ = 0, auv + avwx + u^ + uwx = 0, 

buvz + cuy + cf ^ — cvz = 0, buxz — bcx + cuy + cx^ = 0, 

—b'^u — bcvx + bux — cuy = 0, —cuy — cvxz + uvz — uz^ = 0, 

bwxz + cwy + cxz + cz ~ c = 0, abvz + acy + bcv + be — c = 0, 

—auy — avxz — uxz ~ ux + u = 0, —buv — bvwx — uv — uwy + u = 0, 

—bcu + bu'^z — c^vx + CUV + cux — cuz = 0. 

By just considering the first four equations the problem can be split into 9 different cases, 
and each one of them can then be solved rather easily. After eliminating those solutions 
that reduce to 2-simplex solutions and those with noninvertible A we find 3 basic solutions 
from which others are obtained by the allowed transformations. These solutions and their 
nonhomogeneous additions will be discussed below. 

10 



5.1 



The first base solution is 



A 



(1) 



1 

1 




1 
d 



and when inhomogeneities are added it sphts into two: 



[A\B]^t^ 



1 -1 


X 




1 1 


y 


, [A\Br' = 


1 


_ 





1 

1 




1 
d 



For [A|-B]3 X oi y can be still eliminated by a gauge transformation, for [A|-B]3 the 
gauge freedom has already been used above. 

The transpose of A^ is not obtained by central reflection and therefore constitutes 
another solution: 



[A\B] 



{Ita) 



1 


" 


1 





-1 1 1 


z _ 



, [A\B] 



(itb) 



1 


" 


1 





-d 1 d 


_ 



These forms cannot be changed by gauge, except hj z -^ uz. 



5.2 



There are two upper triangular solutions 



A 



(2) 



a 1 — ah a{bc — 1] 
b I -be 

c 



A 



(2t) 



a 1 — ab c{ba — 1] 
b I -be 

c 



They differ only in the upper right hand entry and are related by transposition and 
central reflection (followed by a <-* c) . But since transposition is not a symmetry of the 
inhomogeneous part they have to be analyzed separately. 

Depending on which parameters have unit value we get three solutions: 



[A|i?]f) 



1 1-6 6-1 


X 




6 1-6 


-bz 


, [AB]f^ = 


1 


z 





a 1 — 6 a{b — 1] 
6 1-6 

1 



abz 
—bz 

z 
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[A\Bt^ 



a 1 — ah a{bc — 1) 
b 1-bc 

c 

This solution illustrates nicely how added freedom in A decreases freedom in B. 

(2t) 

For the transpose A^ we get two solutions 

a 1 — ab c{ba — 1) 

b 1-bc 

^00 c 

5.3 

For the next solution the inhomogeneous terms can always be gauged away and we have 





[Ali?]f ") = 


a I — a a — 1 
1 


X 

y 


, [AB]?'^ = 




1 


_ 





[A\B 



(3) 



a 
1-ab b 1-bc 
c 

Here we have assumed that at least one of a, 6, c is 7^ 1, else we get a diagonal solution 
with arbitrary B. 

The transpose is again a separate case, we get first 

a 1 — ab 
6 
1-bc c 



[A\B] 



(Si) 



-ay 

y 

-cy 



Now if 6 7^ 1 the inhomogeneous part can be eliminated, and we have finally two solutions 



[A\B] 



(3ta) 



a 1 — a 


-ay 


1 


y 


1-c c 


-cy _ 



[A\B] 



(3tb) 



a 1-ab 


' 


b 





1-bc c 






Note how this case is built up from 2-simplex solutions [A|i?]'^^^\ but not as simple tensor 
products. 

5.4 D = 2 

When all indices are modulo 2 only two solutions remain (in addition to reducible ones) 
namely 



"oil 


" 




"010 


" 


1 1 


y 


5 


1 





1 


_ 




1 1 1 


y . 
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where ?/ G Z2 and we have used the gauge freedom to ehminate Bi in the first case. Here 
it might be useful to record the corresponding i?-matrices for y = 0: 



fl . . . 




. . 1 . 




. 1 . . 




... 1 






... 1 




. 1 . . 




. . 1 . 


\ • • • • 


1 . . . ^ 



/ 



These bear some resemblance with know solutions 



\ 



6 Results for the 4-simplex equation 

For the 4-simplex case the 4x4 index matrix is embedded into 10 x 10 matrices in 
four ways. The equations resulting from (^ were solved using the groebner-package 
of REDUCE [|T3|. From the results we eliminated those solutions for which A was in 

Xajl 1\/TJ2J3J4 



?J1J2J3J4 



a block form corresponding to i?'s with tensor products form Rili^t^^^ 



pjlj2j3j4 

-"'«ii2J3M 



]\/TJi32J3 SjO-H T}n3233H 



Ty-jlJ2TJ3J4 



where M is a solution of the 3-simplex 
equation and K, L of the 2-simplex equation. From the remaining list we eliminated all 
cases obtained from the basic ones by central reflection or by inverse, and those with 
singular A. Furthermore we considered only the generic case of a prime D. 

The solutions A of the homogeneous equation (and their transposes) were next used as 
starting points for constructing the non-homogeneous part B. Then the continuous gauge 
freedom was applied to eliminate some freedom from B. The final result is as follows: 



Permutation blocks 



6.1 






1 





-1 


bi' 


1 





1 





62 











1 











1 





_ 
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After a gauge transformation (^TD we would get 5* = [bi + v,b2 — v,0,0] and we could 
eliminate either 61 or 62- 



6.2 

^010 

10 

—a 1 a 1 — ab 

b 

In the generic case we get -B* = [0, 0, z{b — 1), —z{ab — 1)] but this can be eliminated by 
the gauge transformation. Only if a = 6 = 1 would we get something that cannot be 
gauged away, but in that case the system reduces to a 3-simplex solution. 



6.3 



The transpose of the above solutions is a separate case, and yields 



1-10 
10 10 
10 
1-6 b 



bi 
62 





1 -a 

10 10 

a 

1-ab b 



(where 61 or 62 could be gauged away). Again when the A part is restricted the B part 
gains some freedom. 



6.4 



The next cases are somewhat similar to the above, we get 



6.5 



-a 


a-1 


61 ■ 




1 





62 




a 


1-a 





; 





1 


_ 





1 

1 





10 

10 

—a 1 a 

a-1 1-a 1 



-a 


ab-1 


' 


1 








a 


1-ab 








b 


_ 



1 

1 





10 

10 

—a 1 a 

a6 — 1 1 — a6 b 
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6.6 



1 

1 1 
0a 



a a 
-1 
1 — a 







ax 

—X 

—ax 

X 



1 —a ah 

10 1 -h 

a 1 — ah 

h 



6.7 



For the transpose of the above nothing 


can be gauged away, and we g 


et 




" 10 







" 1 









1 







10 









-1110 


X 


1 


-1110 


—hx 


1 




1-10 1 


y 




h -6 1-6 6 


X 




r 


10 


' 




" 1 





-1 




10 







10 









—a I a 





' 


-a I a 









a —1 1 — a 


1 


X 




ah —6 1 — ah 


6 








6.8 

The next A matrix is invariant under central reflection, and gauge transformation changes 
nothing. We get three different cases 



1 1 

1 
-1 



-1 





X 




1 



















; 


1 


1 


y . 





1 











-1 


" 




) 


1 















•) 


d 


1 d 


y . 





a 1 



1 




d 



-a 





" 


1 

















1 


d 


_ 



6.9 



For the transpose of the above the inhomogeneous part is quite different and we get 



1 











1 





1 


— 


-1 


1 





1 











1 



X 

y 

—X 

-y 



a 











1 





1 


— c 


-a 


1 





1 











d 



In both cases there are two free parameters. 
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6.10 



a 


i 


—a 


a 








1 


— 





1 





1 











1 



a 


1 


—a 


ad 


' 








1 


-d 








1 





1 














d 






Triangular blocks 



6.11 



1 l-h h-l l-h 

h l-h h-l 

10 

1 



6.12 



1 

l-h h 

h-l l-h I 

l-h h-l 1 



X 

-hx 

y 

z 



I 











l-h 


h 








c{h-l) 


I -he 


c 






c{l-h) {he -I) l-c I 



X 

—xh 
chx 

y 



6.13 



a I — ah c{ah — 1) cd{l — ah) 

h I -he d{hc - 1) 

c l-cd 

d 



I 

l-h h 

h-l l-h I 



d{l-h) d{h 



a 

I -ah 

c{ah — I] 



d d 



X 

-hx 

y 

-dy 





h 
I -he 



cd{l — ah) d{hc — I) 









" 











c 








l-cd 


d 


_ 



The detailed analysis of this case leads to some subcases that are identical to those of 
central reflected (|6.12|) and are not repeated here. 



a I — ah a{h — 1) a(l — h) 

h l-h (6-1) 

1 

1 



a I — ah a{hc — I) ad{l — he) 

h I -he d{he - 1) 

c l-ed 

d 
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6.14 



6.15 



6.16 



6.17 



6.18 



a 1 — ab 

6 

1 — be c 1 — cd 

d 



a 

1-ab b 

c{ab — 1) 1 — be c 1 — cd 

d 



a 1 — ab a{b — 1) 

b 1-6 

10 

1-d d 



' 1 


1-6 


6 


-1 





X 




a 


1- 


a 


a-1 








6 


1-6 





-bz 







1 


y 








1 





z 


•) 





1 











1-d 


d 


-dz _ 







1-d d 











a 


1-ab c{ab-l) 





" 













6 1 -6c 



















c 



















( 


) 


1- 


cd 


d 










abx 
—bx 

X 

—dx 



a 1 — ab aibc — 1) 

6 l-6c 

c 

1-cd d 



1 

1-6 6 

be — 1 1 — be c 1 — c 

1 



X 

—xb 
-x{l — be) — cy 

y 



a 

l-a6 6 

a(6c— 1) 1 — be e 1 — ed 

d 
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6.19 



Here and in the following case we have rational entries in the index matrix. 



a 1 — a 

1 

10 

d-1 1-d d 



a 1 — ah 

6 

1/6 

d-h 1-d/b d 



6.20 



a 

1-ab b d-b 

1/b {b-d)/b 

d 



6.21 D = 2 

For D = 2 we have the following new solutions 



10 1 


' 




10 10 







1 
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10 


_ 
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1 


1 


" 


1 





1 


1 











1 

















1 


_ 



1110 


" 




10 







10 





; 


111 


_ 





1 


1 


1 


1 


" 








1 


1 








1 





1 














1 


_ 



7 Discussion 

In presenting this (complete) set of linear permutation type solutions we hope that some of 
them could be used in other studies. These applications may require further conditions, 
but we believe that permutation type solutions are so benign that they should satisfy 
these conditions, if just independence of any spectral parameter is acceptable. 

Another hope is that the solutions can teach us something about the equations them- 
selves. Oner observation in that direction is that some of the solutions fall into patterns 
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{2) ('it) I I 

that seem to continue to any n. For example, A2 , A^ and |6.14| start a pattern that 
seems to continue as 

fli 1 — aia2 . . 

aa .. 

1 — 02(23 03 1 — 0304 

04 

1 — 0405 as 



This band structure could make sense even as an infinite matrix, and perhaps we should 
soon start to think what kind of object the "00-simplex" equation might be. 
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